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Introduction 

Les  caract^ristiques  des  dcoulements  avec  decq£fv__ 
lenient  ont,  depuis  des  ann6es,.fait  l’objet  de  rechejjiJtJ 
ches  analytiques  et  exp^rimentales.  Bien  qu’impor- 
tantes  du  point  de  vue  de  leurs  applications  aussi 
nombreuses  qu’essentielles,  il  est  etonnant  de  cons- 
tater  que  les  connaissances  acquises  sur  de  telles 
zones  de  d^collement  en  soient  encore  4  un  slade 
quasi  £l£mentaire. 

Une  zone  de  d^collemenl  peut  Stre  decrite  en 
consid4rant  un  dcoulement  primaire  et  un  £coule- 
nient  secondaire  ou  induit,  chacun  dtant'  intime- 
ment  li£  k  l’autre  par  l’interm<5diaire  d’une  zone 
dans  laquelle  les  tensions  tangentielles  sontdlevees. 

Le  rile  joud  par  une  telle  zone  de  transition  est 
d’un  intdrdt  fondamental,  car  c’est  k  I’intdricur 
mdme  de  cette  zone  qu’a  lieu  la  partie  ia  plus  im- 
portante  de  la  production  de  turbulence,  caracle- 
ristique  de  tels  dcoulements. 

Cos  quinze  dernieres  anndes  ont  connu  un  tra¬ 
vail  de  recherche  considdrable  sur  les  zones  de 
ddcollcmcnt,  travail  entrepris  sous  la  direction  du 
prot'esseur  H.  Rouse  4  l'lnstitut  de  Rechcrches 
Hvdrauliques  d’lowa. 

Debutant  avec  la  these  de  doctoral  de  H.  C.  Hsu 
[1]  sur  un  dlargissement  brusque  de  conduilc  a 
deux  dimensions,  plusieurs  dtudes  sc  sont  succedd. 

La  premiere  d’entre  elles,  fut  une  these  de  «  Master 
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dwHWieristics  of  flow  with  separation  have 
tj\e  subject  of  analytic  and  experimental  re- 
DDC'igtfrch  for  years.  It  is  one  of  the  many  problems 
about  which  little  is  vet  known,  yet  about  which 
much  should  be  learned  because  of  the  numerous 
as  well  as  important  applications  of  separated  flow 
in  engineering  practice.  A  zone  of  separation  can 
be  well  depicted  by  considering  a  primary  stream 
and  an  induced  flow  intimately  hound  together  by 
an  intermediate  zone  of  high  intensity  of  shear. 
The  role  played  by  such  a  transition  zone  is  of 
primary  interest,  for  il  is  within  this  zone  that 
most  of  the  turbulence  production  characteristic  of 
separated  flows  takes  place. 


These  last  fifteen  years  have  known  a  conside¬ 
rable  amount  of  work  undertaken  under  the  guid¬ 
ance  of  Professor  H.  Rouse  at  lhc  Iowa  Institute 
of  Hydraulic  Research.  Beginning  with  a  doctoral 
thesis  by  H.  C.  Hsu  [1J  concerned  with  a  two- 
dimensional  abrupt  expansion,  several  investiga¬ 
tions  have  followed  in  sequence.  The  first  of  them 
was  a  master’s  thesis  presented  by  M.  Aric  [2] 
about  the  flow  past  a  plate  with  a  continuous 
boundary;  the  second  and  the  third  were  two  llicscs 
concerned  with  the  turbulence  characteristics  of 
the  hydraulic  jump,  presented  bv  S.  Nagaratnam 
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of  Science  »  prdsent£e  par  M.  Arie  [2]  sur  T^coule- 
ment  autour  d’une  plaque  perpendiculalre  en  con- 
iuct  avec  une  paroi  continue.  La  seconde  et  la  troi- 
siimie  furent  deux  th6ses  de  <  Master  of  Science  » 
sur  les  caractdristiques  de  la  turbulence  dans  le 
ressuut  hydraulique,  prdsentdes  par  S.  Nagaratnani 
[3]  et  T.  T.  Siao  [4],  rassembldes  et  discuses  dans 
lu  rdfdrcnce  [5].  La  quatridme  fut  une  thdse  de  doc- 
torat  par  M.  C.  Chaturvedi  [8]  sur  des  dlargisse- 
mcnls  de  conduites  avec  diffdrentes  valeurs  de  Tan¬ 
gle  d’dlurgissement.  Gnfln,  pcu  de  temps  aprds  cctte 
dernidre,  ce  fut  une  thdse  de  doctorat  de  T.  Car- 
mody  [7]  sur  l’Otude  d’un  disque  immergd. 

Le  but  essenticl  de  toutes  ces  Etudes  a  did  de  me- 
surer  les  caractdristiques  moyennes  aussi  bien  que 
turbulentes  et  d’en  dtablir  une  correlation  ensuite, 
au  moyen  des  dquations  de  la  quantltd  de  snouve- 
nient  et  de  Pdnergie.  Par  les  dquations  de  la  quan- 
titd  de  mouvement  et  de  Tdnergle  citdcs  ci-dessus, 
les  caractdristiques  du  champ  de  forces  et  de  Tdcou- 
lement  ont  dtd  dvaludes.  Cependant,  par  la  mdthode 
d’intdgration  utilisde  pour  les  diffdrentes  sections  le 
long  de  Tdcoulement,  tous  les  ddtails  internes  dis- 
paraissent.  Une  analyse  «  ligne  par  ligne  >  plutftt 
que  «  section  par  section  »  semblalt  alors  logique. 

A  la  lumidre  des  diffdrentes  publications  de  Rou¬ 
se  [8,  9,  10],  les  termes  du  thdordme  de  Bernoulli 
pour  un  dcoulement  turbulent  ont  maintenant  dtd 
dvaluds  dans  la  prdsente  dtude  pour  un  dcoulement 
normal  k  une  paroi,  pour  des  dlnrgissements  de 
conduite  et  pour  Tdcoulement  autour  d’un  disque. 
Les  raisons  de  ce  choix  sont  doubles.  D’abord,  loin 
(Tdtrc  diffdrents,  les  dcoulements  dans  des  dlargis- 
sements  de  conduite  sont  en  beaucoup  de  points 
semblables  fi  ceux  que  Ton  rencontre  dans  l’dtude 
des  corps  immergds;  tandis  que  les  pertes  de  charge 
sont  considerdes  pour  le  premier,  su  contrepartie, 
la  rdsistance,  est  considerde  pour  le  dernier,  tous 
les  deux  dtunt  ia  consdquence  directe  du  tourbil- 
!<m.  Ensuite,  les  mesures  ddtaillees  des  caractdris¬ 
tiques  de  Tdcoulement  moyen  ainsi  que  celles  de 
Tdcoulement  turbulent  dtaient  disponibles  pour  une 
telle  analyse  dans  les  dtudes  de  Arie,  Chaturvedi 
ct  Carmody. 

L’dtude  analytique  prdsentde  ici  est  une  tentative 
pour  obtentr  les  caractdristiques  ddtailldes  du  tour- 
billon  duns  chacun  des  cas.  Plus  que  d’etre  une 
simple  vdrifleation  d’une  analyse  globale,  le  thdo- 
rdme  de  Bernoulli  pour  un  dcoulement  turbulent, 
consider^  comme  un  moyen  de  vdrifler  les  mesures 
obtenues,  sera  sous  certaines  hypotheses  un  moyen 
appreciable  pour  la  ddterminatlon  de  la  tension 
tangenticlle  de  turbulence  en  grandeur  ct  position. 


Analyse  *fu  probltma 

ASI'KCT  PHYSIQUE  : 

Le  ddcollement  duns  un  dcoulement  flulde  se  pro- 
dull  soit  Ic  long  d’une  paroi  rectillgne  ou  progres- 
sivemenl  courbde  en  prdsence  d'un  gradient  de 
pression  contrnirc,  soit  ft  un  changcment  brusque 
de  gdomdtric  de  la  paroi.  L’dcoulement  rdsultunt 
du  ddcollement  est  pur  nature  non  uniforme.  Si, 
de  plus,  la  gdometrie  des  parols  restunt  inchungde, 
le  iiombre  de  Reynolds  considdrd  est  tel  que  Tdcou- 
lemont  devienne  turbulent,  il  sera  aussi  non  per- 
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[3]  and  T.  T.  Siao  [4;,  and  discussed  and  gathered 
in  Reference  [5]  ;  the  fourth  wus  n  doctoral  study 
carried  out  by  M.  C.  Chaturvedi  [8]  on  pipe  tran¬ 
sitions  with  different  valu ’s  of  the  expansion  angle, 
this  iatter  being  foiiowed  shortly  after  by  the  doc¬ 
toral  dissertation  by  T.  Curmotly  [7]  on  a  tiisk 
fully  surrounded  by  fluid. 

All  of  these  projects  have  been  much  concerned 
with  measuring  the  characteristics  of  the  mean 
flow  and  of  the  turbulence  and  then  correlating 
them  in  accordance  with  the  momentum  and  energy 
equations.  Through  the  aforementioned  moment¬ 
um  and  energy  equations,  the  characteristics  of  the 
force  field  and  of  the  flow  huve  been  evaluated. 
However,  bv  the  nielhod  of  integration  over  sec¬ 
tions,  all  the  internal  details  were  necessarily  lost. 
A  “line  by  line”  analysis  'bus  seemed  necessary. 

In  the  light  of  Rouse’s  papers  [8,  9,  10],  the 
pertinent  terms  of  the  Bernoulli  theorem  for  tur¬ 
bulent  flow  have  now  been  evaluated  line  by  line 
for  flow  past  a  normal  wall,  for  flow  in  pipes 
with  90°  and  15°  expansions,  ar.d  for  flow  around 
a  circular  disk.  The  reasons  for  the  above  choices 
are  twofold.  First,  far  from  being  different,  flow 
in  a  pipe  expansion  is  very  much  the  same  us  flow 
past  an  immersed  body  in  that,  whereas  the  head 
loss  is  considered  for  the  former,  its  counterpart 
the  drag  is  considered  for  the  latter,  both  being  the 
direct  consequence  of  a  standing  eddy.  Second,  the 
detailed  data  of  the  mean  flow  together  with  those 
of  the  turbulence  were  readily  at  hand  for  such  un 
analysis  in  the  work  of  Arie,  Chaturvedi.  und  Cur- 
mody. 

The  present  analytical  work  is  an  ultciupt  to 
obtain  the  detailed  characteristics  of  the  stable  eddy 
in  each  case.  More  than  a  mere  check  of  the  over¬ 
all  analysis,  the  Bernoulli  theorem  for  turbulent 
flow,  as  a  check  on  possible  local  discrepancy,  will 
be  under  certain  assumptions  a  valuable  tool  in  the 
determination  of  both  the  magnitude  and  distribu¬ 
tion  of  the  turbulent  shear. 


Analysis  of  ths  problam 

Physical  aspects: 

Separation  in  fluid  flow  is  expected  to  occur 
either  along  a  straight  or  smoothly  curved  bound¬ 
ary  in  the  presence  of  an  adverse  pre assure  gra¬ 
dient  and  boundary  shear,  or  at  un  ubrupt  change 
in  boundary  geometry.  Flow  resulting  from  sepa¬ 
ration  is  nonuniform  in  nature.  If.  moreover,  the 
boundury  geometry  being  unchanged,  the  Reynolds 
number  under  consideration  is  such  tliut  tuibu- 
lence  will  occur,  the  flow  will  ulso  be  unsteady,  ui 
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manent,  du  moins  pour  autant  que  l’on  considire 
le  mouvement  turbulent. 

La  prdsente  dtude  illustre  un  cas  typique  de 
ddcollement  apparaissant  4  un  changement  brusque 
de  la  g^omEtrie  de  la  paroi  (k  1’arEte  d’un  disque 
et  derriire  une  paroi  perpendiculaire  &  l’Ecoulement 
1  et  au  ddbut  d’un  Elargissement  brusque  de  con- 

duite).  Sans  d^collement,  la  ligne  de  courant  au 
voisinage  de  telles  aretes  devrait  avoir  une  cour- 
bure  de  valeur  inflnie,  ce  qui  entrainerait  uue  va- 
leur  inflnie  du  gradient  de  pression.  Cette  situation 
irrdelle  nous  conduit  k  la  conclusion  que  le  ddcol- 
lement  doit  se  produire.  En  fait,  si  le  rayon  de 
courbure  de  la' ligne  de  courant  de  separation  est 
petit,  comme  dans  le  cas  du  disque  ou  de  la  paroi 
verticale,  le  point  oil  le  ddcollement  se  produit  sera 
un  point  oil  rdgnera  une  pression  faible.  Tout  dcou- 
lement  avec  ddcollement,  au  moins  dans  les  cas  k 
deux  ou  trois  dimensions,  est  dtroitement  lid  k  la 
formation  d’une  zone  de  courant  de  retour  conte- 
nue  par  le  tourbillon. 

A  l’intdrieur  du  tourbillon  et  de  la  zone  de  dif¬ 
fusion,  l’dcoulement  est  essentieliement  rotationnel 
mais,  malbeureusement  pour  l’analyse  mathdmati- 
que,  avec  une  distribution  inconnue  de  la  vorti- 
cite  [9],  En  dehors  de  la  zone  de  mdlange,  dii  k  la 
turbulence,  la  valeur  de  la  somme  des  energies  cine- 
tique  et  de  pression  est  approximativement  cons- 
tante,  montrant  ainsi  des  caractdristiques  d’irrota- 
tionalitd  (l’bcoulement  dtant  prdsumd  dtre  irrotatio- 
.  nel  en  amont  du  point  de  ddcollement).  La  zone  de 

separation  du  tourbillon,  d’autre  part,  est  caraetdri- 
sde  par  un  gradient  de  vitesse  dlevd  dans  la  direc- 
t  tion  radiale,  qui  k  son  tour  provoque  une  intensity 

dlevde  de  tension  tangentielle  turbulente.  A  l’aval, 
cette  tension  joue  le  double  rdle  d'accdldrer  le  fluide 
dans  le  tourbillon  et  simuitandment  de  le  ralentir 
dans  le  courant  principal.  A  l’aval  du  niveau  du 
tourbillon,  par  trausfert  d’dnergie,  cette  tension 
rend  l’Ecouleroent  plus  uniforme,  rdduisant  la  vi¬ 
tesse  dans  la  partie  supdrieure  et  l’augmentant  dans 
la  partie  infdrieure. 

Formulation  et  interpretation  des  Equations 

FONDAMENTALES  1 

L’auteur  a  utilisd,  dans  les  paragraphes  suivants, 
les  Equations  ddrivdes  ainsi  que  les  termes  employes 
dans  sa  rdfdrence  [8],  En  prenant  comme  base  le 
systdme  de  coordonndes  cylindriques  ddflni  par  i, 
r  et  0,  dans  les  directions  axiaie,  radiale  et  angu- 
laire  respectivement,  soil  u,  v,  w  les  composantes 
respectives  de  la  vitesse. 

Pour  le  cas  du  disque  et  de  l’dlargissement  de 
,  conduite,  les  dquations  de  mouvement  et  de  conti¬ 

nuity  peuvent  dtre  simplifies  par  le  fait  que  l’dcou- 
lement  moyen  est  permanent  d'une  part  et  que 
rdcoulement  est  axisymdtrique  d’autre  part  (la 
composante  tangentielle  de  la  vitesse  moyenne,  de 
meme  que  les  ddrivdes  partielles  par  rapport  k  0 
des  quantities  moyennes  et  des  quantitds  dont  on  a 
pris  la  moyenne  sont  identiquement  nulles). 

Dans  ces  conditions,  et  en  sdparant  les  compo¬ 
santes  de  la  vitesse  en  leurs  parties  moyennes  et 
de  fluctuation  : 

(U  =  Tl  -|-  It',  u  ~  u  +  v’, 


least  as  far  as  the  turbulent  motion  is  concerned. 

The  present  thesis  illustrates  a  typical  case  of 
separation  at  an  abrupt  change  in  boundary  geo¬ 
metry  (i.e.,  at  the  edge  of  the  circular  disk,  at  the 
edge  of  the  wall,  and  at  the  beginning  of  the  abrupt 
expansion).  Without  separation  the  streamline  at 
the  edge  should  have  an  infinite  curvature  accomp¬ 
anied  by  an  infinite  pressure  gradient.  This  unreal¬ 
istic  situation  leads  to  the  conclusion  that  separa¬ 
tion  must  occur.  As  a  matter  of  fact,  if  the  radius 
of  curvature  of  the  limiting  streamline  is  small,  as 
in  the  case  of  the  disk  or  of  the  wall,  the  initial 
separation  point  will  possess  a  low  pressure. 

Flow  with  separation,  in  two-dimensional  and 
axisymmetric  cases  at  least,  is  closely  associated 
with  the  formation  of  a  zone  of  reverse  flow  em¬ 
bodying  the  standing  eddy.  Within  the  eddy  and 
the  diffusion  zone,  the  flow  is  essentially  rotational 
but,  unfortunately  for  a  mathematical  analysis, 
with  an  unknown  distribution  of  vorticity  [9]. 
Outside  the  zone  of  turbulent  mixing,  the  flow  is 
seen  to  retain  approximately  a  constant  sum  of 
pressure  and  velocity  head,  thus  showing  irrota- 
fional  characteristics  (the  oncoming  flow  presumab¬ 
ly  being  irrotational).  The  border  of  the  eddy,  on 
the  other  hand,  is  characterized  by  a  high  velocity 
gradient  in  the  radial  direction,  which  in  turn 
causes  a  high  intensity  of  shear.  Farther  down¬ 
stream,  this  shear  plays  over  the  whole  pattern  of 
the  flow  the  twofold  role  of  accelerating  the  fluid 
in  the  eddy  and  simultaneously  decelerating  that 
in  the  main  stream.  Past  the  zone  of  the  eddy, 
through  transfer  of  energy,  the  shear  will  be  res¬ 
ponsible  for  rendering  the  flow  more  uniform,  red¬ 
ucing  the  velocity  in  the  upper  part  and  increas¬ 
ing  it  in  the  lower  part. 


Formulation  and  interpretation  of  thf.  basic 
equations: 

The  writer  has  made  intensive  use  in  the  follow¬ 
ing  paragraphs  of  the  equations  derived  and  of  the 
technical  terms  used  in  Reference  [8]. 

With  reference  to  the  cylindrical  coordinate  sys¬ 
tem  defined  by  x,  r,  and  0  in  the  axial,  radial,  ang¬ 
ular  directions,  respectively,  let  u,  v,  and  w  be  the 
corresponding  components  of  the  velocity.  For  the 
cases  of  the  circular  disk  and  of  the  pipe  expansion, 
the  equations  of  motion  and  of  continuity  can  be 
simplified  by  the  fact  that  the  mean  flow  is  steady 
on  the  one  hand,  that  the  flow  is  axisymmetric- 
on  the  other  hand  (i.e.,  the  mean  tangentiel  velo¬ 
city  component,  together  will  all  the  partial  deri 
vatives  with  respect  to  6  of  all  averaged  quantities, 
is  identically  nil). 

Under  these  conditions,  and  by  separating  the 
velocity  components  into  their  usual  mean  and  flue 
tuating  parts  (i.e.,  u  =  B  -f  u',  v  =  u  +  v\  w  =  u 


w  =  w  +  u>'). 
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!es  Equations  de  continuity  pour  le  mouvement 
raoven  et  turbulent  s'Ecrivent  : 


-+-  the  continuity  equutions  for  mean  and  tur¬ 
bulent  motion  can  be  expressed  us: 


3  ( ru )  ,  3  (rv) 

~3x  1  dr  =° 

a  (ru'j  a  (r V)  a  (rw')  _  n 


tandis  que  les  Equations  de  mouvement  s’Ecrivent 
[8,  eq.  (14)]  : 


whereas  the  equations  of  motion  assume  the  follow¬ 
ing  form  [8,  Eq.  (14)  ] : 


3a  .  /  du  ,  _  du  .  3 u  \  3  t  d2u 

p  IT  +  p  ( u  aF  +  y  "aT  +  w  Tsf)  ~~  zx  ^  +  pn>  +  »  (i*  + 


32u  ,  _1_  du  32n 
dr*  ‘r  r  3r  +  /-^e- 


3a  .  /  3  a  .  3a  .  3u  in2 

"  31  +P(“  3?+"  37+a,  739  ' 


/  3m 

p  (u  ~aF 


_,^<p+p0)+„(^  +  ^  +  jr^_^  + 


32o _ 2^  div_ 

r2 ae2  r  rae 


3  ^  f d2w  ,  a2u>  ,  1  3m  m  ,  32w 

~  TsF  (P  +  W  +  +  -0FT  +  T'ar'-TT  + 


4.  3.*!L\ 
+  r  r36  / 


Si  la  pression  dEcomposEe  en  valeur  moyenne  et 
de  fluctuation  (/>  =  p  + />'),  et  si  les  diffE- 
rentes  composantes  instantanEes  de  la  vitesse  sont 
introduites  dans  liquation  [3],  aprEs  avoir  pris  les 
moyenues  temporelles,  nous  obtiendrons  les  Equa¬ 
tions  de  Reynolds  [8,  Eq.  (16)]. 


If  the  pressure,  expanded  in  terms  of  its  mean 
and  fluctuating  parts  (i.e.,  p  ~p  -f  p'),  as  well  as 
the  different  instantaneous  component  of  the  velo¬ 
city,  is  introduced  into  Equation  (3),  after  averag¬ 
ing  with  respect  to  time  the  Reynolds  equations 
will  be  obtained  [8,  Eq.  (16)] : 


t  -  du 


- : 0  +*  if + ^  -ir  + -  -is  -  ~  -t w + *» + >■  (-g-  +  -P-  +  tIt 


( —  an  j  —  3a  .  ,  at/ 

pv“  a ?+■’ aT  +“  ■sr+1' 


^  r30  r  . 


- gJT  (?  +  Pfl)  +  g. 


a2?  .  d3v 
dx 2  +  3r2 


,  1  du  _  _v 

"P"  r  3r  ‘  r2 


Nous  pouvons  tirer  inunEdiatement  de  l’Equation 
(4)  1’equation  d’Energie  pour  le  mouvement 
moyen;  aprEs  multiplication  de  chacune  des  Equa¬ 
tions  (4)  et  aprEs  addition  de  ces  trois  Equations, 
nous  obtiendrons  :  [8,  Eq.  (26)]. 


Readily  obtainable  from  Equation  (4)  is  the  work- 
energy  equation  for  the  mean  motion;  multiplica¬ 
tion  of  each  of  Equation  (4)  by  the  corresponding 
mean  velocity  component  and  then  addition  of  the 
three  equations  will  give  [8,  Eq.  (26)]  : 
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La  signification  dc  tous  les  termes  contenus  dans 
cctte  equation  cst  d’un  interet  primordial,  car  de 
cctte  manierc  une  vue  plus  approfondie  du  mEca- 
nisme  fondainental  du  mouveinent  fluide  peut  Etre 
j  obtenue.  Cette  Equation  a  ete  arrangee  de  facon 

I  &  montrer  dans  les  termes  de  gauche  la  variation 

7  du  travail  fait  et  dans  ceux  de  droite  la  variation 

correspondanle  de  I’energie  cinetique  moyenne.  Le 
?  premier  terme  dans  le  membrc  de  gauche  repre- 

|  sente  le  travail  total  fait  par  les  forces  de  press'  in 

j  et  les  forces  de  pesantcur;  le  second,  le  travail  cun- 

1  servatif  fait  par  les  tensions  visqueuses;  le  troi- 

|  siEme,  le  travail  total  fait  par  les  tensions  de  Rey- 

‘  nolds;  le  qu&trieme,  le  travail  fait  de  manierc  dis- 

sipatrice  par  les  tensions  visqueuses,  et  le  cinquieme, 
le  travail  fait  par  les  tensions  de  Reynolds  en  pro- 
duisant  la  turbulence.  Par  suite,  ce  terme  represente 
une  perte  pour  le  mouvement  moyen  en  meme 
temps  qu’un  gain  pour  le  mouvement  turbulent. 
C’est  ensuite  sculcmcnt  cn  examinant  sa  contre- 
partie  dans  l’equation  d’energie  pour  la  turbulence 
qu’on  peut  comprendre  plus  en  dEtail  les  diffdrentes 
i  etapes  du  mEcanisme  de  la  turbulence. 

5  L’Equation  d’energie  pour  les  fluctuations  turbu- 
\  lentes  seules  peut  Etre  etablie  en  multipliant  cha- 
cune  des  Equations  de  Navier  et  Stokes  (Eq.  (3)) 
j  par  les  composantes  respectives  de  la  vitesse,  sEpa- 
rEes  en  parties  moyennes  et  de  fluctuations.  En- 
|  suite,  aprEs  avoir  pris  la  moyenne  temporelle  des 
diffErents  termes,  et  aprEs  avoir  soustrait  l’equa- 
tion  d’energie  (eq.  (5)),  nous  avons  ;  [8,  eq.  (27)]. 


The  significance  of  all  the  terms  embodied  in 
this  equation  is  of  primary  interest,  for  in  this  way 
a  deeper  insight  into  the  basic  mechanism  of  fluid 
motion  can  be  gained.  This  equation  has  been  ar¬ 
ranged  in  such  a  way  as  to  show  on  the  left-hand 
side  the  rate  of  work  done  and  on  the  right-hand 
side  the  corresponding  rate  of  change  of  mean 
kinetic  energy.  The  first  term  in  the  member  on 
the  left  represents  the  total  work  done  by  pressure 
and  body  forces;  the  second,  tlie  conservative  work 
done  by  the  viscous  stresses;  the  third,  the  total 
work  done  by  the  Reynolds  stresses;  the  fourth, 
the  work  done,  in  a  dissipative  way,  by  the  viscous 
stresses;  and  the  fifth,  the  work  done  by  the  Rey¬ 
nolds  stresses  in  producing  turbulence.  This  term, 
consequently,  represents  a  loss  from  the  mean  mo¬ 
tion  to  the  turbulent  motion.  Only  by  examining 
its  counterpart  in  the  turbulent -energy  equation  can 
further  understanding  of  the  different  steps  of  the 
turbulence  mechanism  he  gained.  The  energy 
equation  for  turbulence  fluctuations  alone  can  be 
derived  by  multiplying  each  of  the  Navicr-Stokes 
equations— Equation  (3) — bv  the  corresponding 
velocity  components,  expanded  in  their  usual  mean 
and  fluctuating  parts;  then,  after  averaging  has 
been  done  on  the  various  terms,  subtraction  from 
the  foregoing  mean  energy-work  equation — Equa¬ 
tion  (5) — will  give  [8,  Eq.  (27)]: 
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J  On  peut  distinguer  successivement  dans  cette 
|  equation  les  termes  reprEsentant  :  »e  travail  fait 
j»ar  les  fluctuations  de  pression,  le  transfert  d’Ener- 
gie  5  partir  du  mouvement  moyen  et  le  travail  total 
J  fait  par  les  tensions  visqueuses  dans  le  nicmbre  de 
:  gauche,  Dans  lc  nicmbre  de  droite,  on  reconnait  la 
t  variation  de  l’energie  cinetique  de  la  turbulence 
•  dans  les  termes  qui  represented  la  convection  for- 
j  cee  par  le  mouveinent  moyen  et  la  diffusion  par 
|  la  turbulence,  tandis  que  le  dernier  terme  est  la 
|  source  principale  de  perte  de  l’Energie  turbulente 
1  par  dissipation  visqueuse  dans  les  tourbillons.  Ce- 
!  pendant,  puisque  les  diffErents  termes  de  l’Equa- 
j  tion  (6)  ne  sont  pas  encore  mesurables,  cette  Equa¬ 
tion  sera  essentielienient  d’un  interEt  thEoriquc 
dans  cette  Etude. 


Successively,  in  this  equation  one  can  distinguish 
terms  for:  work  done  by  the  pressure  fluctuations, 
transfer  of  energy  from  the  mean  motion,  and  total 
work  by  the  viscous  stresses  on  the  left-band  side. 
On  the  right-hand  side,  one  can  recognize  the  rate 
of  change  of  the  turbulent  kinetic  energy  in  terms 
which  represent  forced  convection  by  the  mean  mo¬ 
tion  and  diffusion  by  the  turbulence,  whereas  the 
last  term  is  the  main  source  of  loss  of  the  turbu¬ 
lent  energy  by  viscous  dissipation  within  the  ed¬ 
dies.  However,  because  the  different  terms  of 
Equation  (6)  are  not  as  yet  measurable,  the  equa¬ 
tion  will  be  only  of  theoretical  interest  in  this 
thesis. 
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LIQUATION  UK  BERNOILLI  EN  tCOL'LEMEST  TURBl - 
LENT  : 

L’equation  de  Bernoulli,  dans  sa  forme  habituelle, 
valable  seulement  dans  le  cas  d’^coulements  non 
visqucux  et  incompressibles,  peut  £tre  etablie  4  par- 
tir  de  l’equation  d'dnergie  pour  l’dcoulement  moyen 
dt)j4  developpee  pour  le  cas  dVcoulenients  visqueux 
et  turbulents. 

En  posant  B'  =  (1/2)  pV2  -{-  p  +  pfi  dans  laquelle 
V  est  la  vitesse  moyenne  totale,  p  la  pression 
moycnne  et  Q.  la  fonction  «  potentiel  dd  4  la  pesan- 
teur  »,  1’dquation  (5)  peut  s’dcrire  sous  la  forme  : 


Bernoulli  equation  for  tfrbi  lknt  flow: 

The  Bernoulli  equation,  which  in  its  usual  form 
is  valid  only  in  the  case  of  mmviscnus  incompres¬ 
sible  flow,  can,  for  convenience,  be  derived  for  the 
case  of  viscous  turbulent  flow  from  the  previously 
stated  work-energy  equation  for  the  mean  flow. 

By  setting  B'  =  )  pY2  -f-  7>  — |-  pfi,  where  V  is  the 
total  mean  velocity;  p  is  the  mean  pressure;  and 
n  is  the  gravitational  potential,  Equation  (ft)  can 
be  written  in  the  following  form: 


.  3B'  .  _  3B' 
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Aprds  multiplication  de  chacun  des  membres  de 
cette  Equation  (7)  par  le  ddplacement  dleinentaire 
le  long  d’une  ligne  de  courant  et  en  modifiant  le 
premier  membre,  on  peut  l’dcrire  sous  la  fonne  : 


After  multiplying  each  side  of  this  Equation  (7) 
by  the  elementary  displacement  ds  along  a  stream- 
tube  and  modifying  the  first  term  through  ruies 
of  calculus,  by  noting  that  u  =  dx/dt  and  T>  =  dr/dt 
on  can  write: 


3B'  ,  „  3B' 
dx 


3B ’\  , 

+"sr)d— 


VdB' 


(8) 


Dans  cctte  Equation  d B'  est  la  diffdrentielle  totale 
de  la  somme  [(1/2)  pV2  +  p  4  pD] .  On  peut  alors 
evaluer  l’intdgrale  de  liquation  (7)  entre  deux  sec¬ 
tions  repdrdes  par  les  indices  1  et  2  pour  donner  : 


Herein  dB'  is  the  total  differential  of  the  sum 
(i  pY2  -f-  p  -f-  pH).  Equation  (7),  now,  can  be  inte¬ 
grated  along  a  streamtube  between  two  sections 
indicated  by  the  subscripts  1  and  2  to  yield: 


Les  dilTdrents  terines  de  cette  Equation  pourraient 
dire  discutds  de  la  nidme  fa$on  que  font  ceux 
de  liquation  (5).  D’une  manure  analogue,  le  terine 
pour  la  production  de  la  turbulence  aura  sa  contre- 
partie  sous  la  forme  d’un  gain  pour  l’dcoulement 
turbulent  aux  ddpens  de  l’dcoulement  moyen. 

Tandis  que  le  thdorcine  de  Bernoulli  dans  la 
fo  me  ci-dessus  est  rigoureuse  en  elle-meiue,  il  est 
cvidein ment  trnp  complexe  pour  Dtre  utilise  pour 

07(5 


The  different  terms  of  this  equation  have  the 
interpretations  already  discussed  for  Equation  (5). 
In  a  similar  fashion,  the  term  for  production  of 
turbulence  will  have  its  counterpart  in  the  form 
of  a  gain  for  the  turbulent  flow  at  the  cost  of  the 
mean  flow. 

While  the  Bernoulli  equation  in  the  above  form 
is  rigorous  n  itself,  it  is  too  complicated  to  he 
readily  used  for  computations.  However,  the  cqua- 
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les  calculs.  On  peut  renurquer  que  liquation 
pourra  etre  grandement  siinplifiee  et  en  fait  rendue 
plus  apte  au  calcul  en  notant  que  duns  les  cas 
,  consider^,  les  forces  dues  k  la  pesanteur  ne  jouent 
aucun  role  et  que  les  tensions  visqueuses  peuvent 
etre  negligees  en  comparaison  des  tensions  de  Rey¬ 
nolds  pour  les  nombres  de  Reynolds  elcves  dont  il 
est  question. 

Liquation  ainsi  simplifiee  peut  s’ecrire  : 


tion  can  be  greatly  simplified  and,  in  fact,  rendered 
more  amenable  to  computation  by  noting  that,  in 
the  cases  under  consideration,  the  gravitational 
forces  play  no  role,  and  the  viscous  stresses  can 
be  neglected  in  comparison  with  the  Reynolds 
stresses  for  the  high  Reynolds  numbers  involved. 
The  simplified  equation  thus  assumes  the  form: 
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Et  liquation  correspondante  pour  le  cas  a  deux  The  corresponding  equation  for  the  two-dimon- 

dimensions  pourrait  &tre  etablie  de  la  meme  ma-  sional  case  could  be  derived  in  the  same  way  to 

nifere  pour  donner  :  yield: 


B'j  ==  B'n  +  p  (it  U'2  V  u’.n  -f  (n  «v  +  Ti  />")  J 
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Discussion  des  rtsultats 

Valeurs  pr£vues  : 

Comine  premier  pas  vers  la  discussion  des  r<5sul- 
tats,  lPquation  de  Bernoulli  sous  forme  adimen- 
sionnelle  doit  fitre  rappePe  des  chapitres  prece¬ 
dents  : 

=  B2  +  T  +  L 

Dans  cette  expression,  il  va  sans  dire,  puisque 
l’energie  est  tou jours  transferee  de  recoupment 
moyen  vers  1’ecoulement  turbulent  et  n’est  jamais 
r6cuperee  par  recoupment  moyen,  que  le  terme  de 
production  de  la  turbulence  L  est  toujours  positif, 
tandis  que  T  ne  conserve  pas  obligatoirement  un 
nPme  signe  puisqu’il  represente  un  transfert 
d’dnergie  d’un  tube  de  courant  &  un  autre. 

Comme  les  autres  termes  de  liquation  de  Ber¬ 
noulli  deji  consider,  le  terme  de  transfert  a  les 
dimensions  d’energie  par  unite  de  volume.  Si,  de 
plus,  la  valeur  du  transfert  est  multipliee  par  le 
volume  par  unite  de  temps  urfA  sPcoulant  k  tra- 
vers  une  section  et  que  l’on  prenne  l’integrale  a 
travers  toute  la  region  de  recoupment,  le  travail 
interne  sur  les  surfaces  des  elements  sPliminerait 
et  il  resterait  seul  le  travail  fait  par  les  tensions 
de  Reynolds  sur  les  sections  extremes. 

Cependant,  ce  fait  ne  nous  fournit  pas  une  veri¬ 
fication  supplemental  puisqu’il  s’agit  simplement 
de  l’equation  d’energie  pour  recoupment  moyen 
dans  sa  forme  integrale  d6j&  utilis^e  et  faite  pr£- 
cedemment  [2,  6,  7]. 

La  sommc  des  trois  termes  B,  L  et  T,  de  toute 
evidence,  doil  posseder  une  valeur  constunte  et 


Discussion  of  results 

Expected  values: 

As  a  first  step  towards  discussing  the  results,  the 
equation  of  Bernoulli  in  its  nondiiuenslonal  sym¬ 
bolic  form  must  be  recalled  from  the  previous 
chapter : 

B,  =  B2  +  T  +  L 

Herein,  needless  to  say,  since  energy  is  always 
transferred  from  the  mean  flow  to  the  turbulent 
flow  and  is  never  recovered  by  the  mean  flow,  the 
turbulence-production  term  L  is  always  positive, 
whereas  T  is  not  restricted  to  one  sign,  since  il 
represents  a  transfer  of  energy  from  one  stream- 
tube  to  another.  Like  the  other  terms  of  the  Ber¬ 
noulli  equation  already  considered,  the  transfer 
term  is  an  expression  of  energy  per  unit  volume. 
If,  moreover,  the  value  of  the  transfer  is  multiplied 
by  the  volume  per  .unit  time  ud A  flowing  across  a 
section  and  integrated  across  the  entire  flow  region, 
the  internal  work  on  the' surfaces  of  the  element 
will  cancel  and  only  the  work  done  by  the  Reynolds 
stresses  on  the  end  sections  could  remain.  How¬ 
ever,  this  fact  does  not  provide  a  supplementary 
check,  since  it  is  simply  the  analysis  of  the  energy 
equation  for  the  mean  flow  in  its  integral  form 
previously  performed  [2,  6,  7]. 

The  sum  of  the  three  terms  B,  L,  and  T,  quite 
obviously,  must  have  a  constant  magnitude  of  unity 
everywhere  outside  the  eddy,  for  the  approaching 
flow  possesses  this  value.  Within  the  eddy,  on  the 
other  hand,  every  streamline  forms  a  loop  and  the 
locus  of  its  origin  for  the  computations  has  been 
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4gale  a  Funite  partnut  en  dehors  du  tourhillon  car 
l’eeoulement  exterieur  voisin  possede  cctlc  valeur. 

A  1’interieur  du  tourhillon,  d’autre  part,  chaque 
ligne  dc  courant  forme  unc  boucle  et  son  origine 
pour  les  calculs  a  etc  arbitraircment  choisie  cornme 
dtant  la  ligne  de  composante  axiale  de  vitesse  nulle. 
Le  long  de  telles  lignes  de  courant  females,  la 
soimne.  B  -f-  L  +  T  sera  par  consequent  constante, 
egule  a  la  valeur  Ho  quo  Ton  a  choisie  a  son  ori¬ 
gine.  La  pression  a  trovers  l’ecoulenient  est  essen- 
tielleinent  egale  4  celle  regnant  en  dehors  du  tour- 
billon  el  le  terme  de  pertes  L  possede  une  valeur 
beaucoup  plus  faible  que  celle  trouvee  dans  1’tScou- 
Icment  principal.  De  plus,  la  sotniuc  des  termes  dc 
Bernoulli,  coinprenant  les  termes  de  transfert  et 
de  pertes  aura  une  valeur  beaucoup  plus  faible  que 
celle  de  l’dcoulement  principal,  &  cause  de  la  valeur 
negative  du  transfert  T. 

II  Jecoule  imnnidiatement  de  ces  considerations 
que  le  terme  de  pertes  L  doit  ndeessairement  etre 
egal  au  terme  de  transfert  T  quand  une  ligne  de 
courant  se  referme  sur  ellc-meme. 

RESULTATS  DES  CALCULS  I 

Mis  a  part  les  cas  d’un  dlargissement  dc  conduite 
de  15°  de  la  figure  4,  du  disque  de  la  figure  6  el  de 
la  paroi  perpcndiculaire  4  l’dcoulement  de  la  figu¬ 
re  8,  les  resultats  sont  en  defaut  dans  les  premie¬ 
res  sections,  tandis  qu’ils  sont  en  exces  dans  les 
sections  suivantes,  la  seule  exception  etant  le  cas 
de  l’eiargissement  brusque  de  90°  presente  dans  la 
figure  2.  Ces  resultats,  si  l’on  fait  l’hypothese  que 
les  caracteristiques  inoyennes  de  l’ecoulement  ont 
dte  correctement  determinees  dans  chaevn  des  cas, 
montrent  que  les  mesures  des  tensions  de  la  turbu¬ 
lence  sont  trop  faibles  dans  les  premidres  sections 
mais  trop  importantes  dans  les  sections  suivantes. 
On  peut  dtre  ainsi  amene  A  penser  que  les  mcsurcs 
dc  la  turbulence  plutol  que  les  mesures  des  vitesses 
et  pressions  inoyennes  ont  etc  incorrectement  dd- 
termindes  [  V  ]  - 

Correction  des  result ats  exp£rimentaux  : 

Si  Ton  fait  l’hypothdse  que  les  caractdristiques 
inoyennes  de  I’dcoulement,  de  meme  que  celles  des 
tensions  normales,  ont  dte  correctement  deterini- 
ndes,  une  premidre  approximation  de  la  valeur  de 
la  tensi>  .!  tangentielle  u'v'  due  4  la  turbulence  est 
obtenue  par  un  calcul  par  approximations  succes- 
sives  en  considdrant  seulement  les  termes  compre- 
nant  liV  multiplid  par  un  gradient  de  vitesse  dans 
une  direction  particulidre.  Une  correction  de  cette 
approximation  est  alors  possible  de  maniere  4  tenir 
compte  de  toules  les  expressions  conlenues  dans 
les  termes  de  pertes  el  de  transfert,  Ce  calcul  a  dte 
fait  dans  cette  etude  pour  le  cas  du  disque. 

Les  valeurs  de  la  tension  tangentielle  de  turbu¬ 
lence  mesurde  prdeddemment  [7]  sont  reprdsentdes 
dans  la  figure  9,  tandis  que  celles  eorrespondant  au 
calcul  sont  representds  sur  la  figure  10.  Ces  valeurs 
ne  sont  pas  ndeessairement  ce  qu’elles  devraient 
dtre,  car  1'hypothdse  que  les  tensions  normales  aient 
dte  correctement  mesurdes  nc  peut  pas  dtre  justi- 
lide. 

En  fait,  il  est  probable  sinon  certain  que  les 
vnleurs  mesurdes  des  tensions  normales  de  la  tur¬ 
bulence  sont  trop  faibles.  En  consequence  de  quoi, 
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chosen  arbitrarily  as  the  line  of  zero  forward  velo¬ 
city.  Along  such  closed  streamlines,  the  sum 
B  4-  L  T  will  thereby  lie  a  constant,  equal  to  the 
value  B„  assumed  at  its  origin.  The  pressure  across 
the  flow  is  essentially  the  same  as  that  outside 
the  eddy  and,  due  mainly  to  the  low  velocity  head 
inside  the  eddy,  the  loss  term  L  is  seen  to  have 
a  much  smaller  vaiue  than  that  of  the  primary 
flow.  Moreover,  due  to  a  negative  value  of  the 
transfer  term  T,  the  value  of  the  Bernoulli  sum, 
including  terms  for  transfer  and  loss,  will  be  seen 
to  retain  a  much  lower  value  than  that  of  tlu*  pri¬ 
mary  flow.  Readily  seen  from  these  considerations 
is  the  fact  that  the  loss  term  L  must  necessarily 
equal  the  transfer  term  T  when  a  streamline  closes 
upon  itself. 


Results  of  the  computations: 

In  the  case  of  the  15°  pipe  expansion  in  figure  4, 
the  circular  disk  of  Figure  6,  and  the  normal  wall 
of  Figure  8,  the  results  are  deficient  in  the  first 
sections,  whereas  they  arc  in  excess  in  the  follow¬ 
ing  sections,  the  only  exception  being  the  case  of 
the  90°  pipe  expansion  presented  in  Figure  2. 
These  results,  if  one  assumes  that  the  mean  char¬ 
acteristics  of  the  flow  have  been  correctly  deter¬ 
mined  in  either  case,  show  that  the  measurements 
^  the  turbulence  stresses  are  too  low  in  the  first 
sections  but  too  high  in  the  next  .  Conse¬ 

quently,  one  may  be  led  to  think  that  the  turbu¬ 
lence  measurements  rather  that  the  mean  pressures 
and  velocities  were  incorrectly  determined  [7]. 


Correction  ok  the  experimental  results: 

If  one  assumes  that  the  mean  characteristics  of 
the  flow,  together  with  tbe  normal  stresses,  have 
been  correctly  determined,  a  first  approximation  of 
the  turbulent  shear  u'v'  is  obtained  from  a  trial- 
and-error  computation  by  considering  only  the 
terms  involving  TiV  times  a  gradient  of  velocity  in 
a  particular  direction.  A  correction  of  this  ap¬ 
proximation  is  then  possible  in  order  to  take  into 
account  all  the  expressions  contained  in  the  loss 
and  transfer  terms.  A  computation  of  this  sort 
has  been  made  for  the  case  of  the  circular  disk. 

The  values  of  the  turbulent  shear  measured 
previously  [7]  are  presented  in  Figure  9,  whereas 
those  required  by  the  analysis  arc  presented  in 
Figure  10.  These  values  are  not  necessarily  cor¬ 
rect,  for  the  assumption  of  correct  determination 
of  normal  stresses  is  not  justified.  Indeed,  it  is 
probable,  although  not  certain,  that  the  measured 
values  of  the  turbulence  normal  stresses  are  too 
low.  Consequently,  the  correct  value  that  the  tur¬ 
bulent  shear  should  assume  will  lie  somewhere 
between  the  curve  of  measured  values  and  that  of 
the  values  required  by  the  analysis. 

This  trial-and-error  computation  shows  that  the 


Ies  valcurs  correctes  que  les  tensions  tangentielles 
dcvraient  avoir  seraient  comprises  entre  les  cour- 
bes  dcs  valcurs  mesuroes  et  de  celles  requises  par 
i’analvse. 

Le  calcul  par  approximations  successives  montre 
quc  le  terme  de  transfert  depend  principalement 
des  valcurs  du  gradient  dans  la  direction  radiale 
du  produit  de  u'v'  par  rii  et  du  gradient  dans  la 
direction  axiale  du  produit  de  T75  par  ru.  Tandis 
que  !e  lerme  de  pertes  depend  plus  des  valeurs 
itiemes  de  u'v'  et  de  li'1.  En  d’autres  termes,  plus  la 
pente  de  la  courbe  representative  du  produit  de  ru 
par  u'v'  en  function  de  r  est  grande,  plus  le  trans¬ 
fert  d’dnergie  est  important. 

En  fait,  on  s’apercoit  que  le  signe  de  ce  transfert 
est  determine  principalement  par  le  signe  de  la 
pente  de  cette  courbe.  Le  terme  de  pertes,  L,  d’autre 
part,  sera  toujours  positif  et  proportionne!  princi- 
palement  aux  valeurs  de  la  tension  tangentielle  de 
turbulence  et  au  gradient  de  la  vitesse  axiale 
moyenne  dans  la  direction  radiale. 


Application  do  liquation 

Application  num£rique  : 

Le  theoreme  de  Bernoulli  derive  dans  le  chapitre 
precedent  n’est  evidemment  pas  pratique  pour  les 
calculs.  D’abord,  l’intcgrafion  le  long  d’une  ligne  de 
courant  est  tres  difficile  a  efTectuer.  Ensuite,  de 
maniere  a  conserver  la  generalite  des  expressions, 
les  equations  (9)  et  (10)  doivent  etre  ecrites  dans 
leurs  formes  adinicnsionnelles. 

En  inlroduisant  la  relation  ds/V  —  dx/u  dans 
l’equation  de  Bernoulli  et  en  divisant  par  la  valeur 
de  l’energie  dans  l’ecoulement  d'approche  pU,,2/2,  il 
en  resulte  la  forme  adimensionnelle  : 
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transfer  term  depends  mainly  on  the  value  of  the 
gradient  in  the  radial  direction  of  the  product  of 
u'v'  by  ru  and  on  the  gradient  in  the  x  direction 
of  the  product  of  7 F5  by  ru,  whereas  the  loss  term 
depends  more  on  the  values  of  the  u'v'  and  u”~ 
themselves.  In  other  words,  the  greater  the  slope 
of  ru  times  7717  versus  r,  the  greatei  the  transfer 
of  energy.  The  sign  of  the  transfer  term  in  indeed 
seen  to  be  determined  mainly  by  the  sign  of  the 
value  of  this  slope.  The  loss  term,  on  the  other 
hand,  will  be  always  positive  and  directly  propor¬ 
tional,  mainly,  to  the  values  of  the  turbulence  shear 
stresses  and  the  gradient  of  the  mean  axial  veto 
city  in  the  radial  direction. 


Application  of  the  equation 

Numerical  procedure: 

The  Bernoulli  equation  derived  in  the  preceding 
chapter  is  evidently  not  convenient  for  computa¬ 
tional  use  for  two  reasons:  first,  integration  along  a 
stream  line  is  difficult  to  handle;  second,  in  order 
to  increase  the  utility  of  the  expressions.  Equa¬ 
tions  (9)  and  (10)  must  be  written  in  non-dimen¬ 
sional  form. 

Upon  introducing  the  relationship  ds/V  —  tlx /ft 
in  the  Bernoulli  equation  and  dividing  by  the  value 
of  the  energy  level  in  the  oncoming  flow,  pl'„'-/2, 
the  following  dimensionless  expression  will  result: 


B,  =  B,,  +  2 
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(12) 


dans  laquelle  B  est  la  contrepartie  sans  dimension 
de  B',  Un  la  vitesse  dans  l’ecoulement  du  courant 
d’approche  et  D„  une  longueur  caracteristique;  dans 
ce  cas,  le  diametre  de  la  section  d’entrde  au  ddbut 
de  Pelargissement. 

D’une  maniere  tres  seinblable,  Tequation  de  Ber¬ 
noulli  pour  un  ecoulement  turbulent  a  deux  dimen¬ 
sions  pent  etre  etablie.  L’equation  adimensionnelle 
simplifide  se  reduit  a  : 


Herein,  B  is  the  dimensionless  counterpart  of  B', 
U0  the  velocity  in  the  approaching  flow  and  D0  a 
characteristic  length,  in  this  case  the  diameter  of 
the  intake  end  of  the  pipe  expansion. 

In  a  very  similar  way,  the  Bernoulli  equation  for 
two-dimensional  turbulent  flow  can  be  derived. 
The  simplified  dimensionless  equation  reduces  to: 


B, 


».  4-  2 
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dans  laquelle  h  est  une  longueur  caracteristique  in  which  h  is  a  characteristic  length  chosen  as  the 

choisic  enmme  etant  la  hauteur  d’une  paroi.  height  of  the  wall. 
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II  i-st  it  noter  que  Ifs  tenues  duns  les  equations 
(12)  ft  (i.ii  sunt  lies  integrales  curvilignes.  Kn  d’au- 
tris  teniifs,  flits  vmt  function  tie  I'abscisse  curvi- 
ligne  */!)„  mi  s  It  If  long  il'uiif  ligne  iif  courant. 
Hn  fiinsi'ij nfiiff,  tnutes  les  quantities  ci-iiessus  out 
ftf  falcu Iffs  dans  la  |iri;si>nlf  etude  pour  loute  la 
region  df  I'l'cmilf nu'iil  on  un  nombre  dc  points 
ausxi  grand  que  possililf,  L<‘s  valours  prises  par  ces 
functions  le  lung  d’unr  ligne  de  oourant  out  etc  tra- 
fffs  fii  function  ties  \aleurs  oorrospondantcs  do 
■r/I)„  im  xh  'la  longueur  do  l’aro  lo  long  d'une  ligne 
de  eotlianl  olant  prise  eoinme  function  des  eoordon- 
noos  r/D„  mi  x/h  pour  I ’inleg ration ). 

hour  tons  le,s  oaleuls,  les  deux  hypotheses  ayant 
trail  a'.x  oarao terisliques  de  ia  turbulence  ct  qui 
hirent  rmisideroos  I'niniiio  justifioes  dans  les  io- 
olierolies  juecodentos  id,  7)  out  eto  utiiisees.  Ces 
hypotheses  stmt  :  o'-  pour  des  angles  de  90  “ 

et  hV  d'elargissoinont  do  conduite  ■  ct 

n'-  —  2  >;*  =-=  2  w"- 

pour  lo  disque. 


It  must  be  noted  that  the  terms  under  the  inte¬ 
gral  sign  in  Equations  (12)  and  (13)  are  line  inte¬ 
grals.  In  other  words,  they  are  functions  of  the  arc 
length  .</Du  or  s/h  along,  a  streamline.  Consequent¬ 
ly,  all  the  above  quantities  have  been  computed  in 
the  present  thesis  throughout  the  flow  region  at  ns 
many  points  as  possible.  The  values  taken  by  these 
functions  along  a  streamline  have  then  been  plott¬ 
ed  against  t lie  corresponding  values  of  .r/I),,  or  x/h 
(the  are  length  along  a  streamline  being  a  func¬ 
tion  of  the  coordinates  ,r/U(1  or  x/h)  for  integration. 

Throughout  the  computations,  the  two  assump¬ 
tions  with  regard  to  the  turbulence  characteristics 
which  were  considered  satisfactory  by  t tic  previous 
investigator  [6,  7]  were  used.  They  arc :  />'-  —  Ti™ 
for  90"  and  la"  pipe  expansion;  u'-  2  o'-  —  2  in'- 

for  the  circular  disk. 


Itepailition  des  tenues  de  Bernoulli,  /2 /  Distribution  of  Bernoulli  terms  for  90-  pipe  expansion 
pom-  un  iuyau  divergent  de  SD" 
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Pit  KS  l-XT  AT  I  OX  UKS  lU-lSl'LTATS  : 

Les  resultats  tie  Pan:! lyse  sont  presentes  sous 
forme  graphiquc.  Les  figures  1,  3,  5  et  7  montrent 
ies  lignes  de  courant  traeees  a  l'exterieur  et  a  l'in- 
terieur  du  tourbillon  et  reperees  par  les  valeurs  de 
la  function  du  courant  4/.  Les  difTerents  tenues  de 
i’equation  de  Bernoulli  sont  represents  dans  (es 
figures  corresjtondanles  2,  4,  6  et  8.  Ils  sont  traces 
necessairement  d’unc  nianiere  cumulative  puisquc 
ce  qui  a  etc  perdu  on  transfere  entre  deux  sections 
est  evidemment  perdu  ou  transfere  pour  les  sec¬ 
tions  suivantes. 

Dans  1’ecoulement  principal,  la  sonnne  adimen- 
sionnclle  de  Bernoulli  possede  une  valeur  egale  4 
l’unitc  ]>our  la  section  d’entree,  cetle  section  a  done 
etc  choisie  coniine  origine  pour  les  calculs.  Cepen- 
dant,  a  1’interieur  du  tourbillon,  aucune  section 
semblable  n’cxiste;  par  consequent,  l’intersection 
de  chaque  ligne  de  courant  avec  la  ligne  lieu  des 
points  de  composante  axiale  de  vitesse  nulle  a  ete 
choisie  arhitrairement  comme  origine  pour  les  cal¬ 
culs. 


Prkskxtatio.v  or  bksilts: 

The  results  of  the  aforementioned  analysis  ;ue 
presented  in  graphical  form.  Figures  1.  ,'i,  a,  and  7 
show  the  patters  of  streamlines  extended  into  tlie 
eddy  zone,  and  referred  to  by  the  values  of  the  usual 
stream  function  ij/.  The  terms  of  tlu  Bernoulli 
equation  are  presented  in  the  corresponding  Fig¬ 
ures  2,  4,  6  and  8.  They  are  necessarily  plotted  in  a 
cumulative  way,  since  what  has  been  losl  or  trans¬ 
ferred  between  two  sections  is  evidently  losl  or 
transferred  for  the  following  sections. 

In  the  primary  flow,  the  non-dimensional  Ber¬ 
noulli  sum  assumes  a  value  of  unity  at  the  intake 
section;  this  section  has  consequently  been  chosen 
as  the  initial  section  for  the  computations.  How¬ 
ever,  within  the  eddy,  no  such  initial  section  exists; 
hence  the  intersection  of  every  streamline  with  t in¬ 
line  of  zero  forward  velocity  has  been  chosen  arbi¬ 
trarily  as  origin  for  the  computations. 
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Dans  eette  etude,  outre  unc  discussion  generate  In  this  thesis,  a  general  discussion  of  the  conver-  t 

lie  la  transformation  de  I’energie  dans  un  ecoule-  sion  of  energy  in  flow  with  separation  has  been  pro*  j 

meal  avec  decol lenient,  la  variation  des  differents  sented.  More  specifically,  the  variation  of  tile  dif-  ; 

tenues  de  l  iquation  de  Bernoulli  le  long  d'une  li-  ferent  terms  of  the  Bernoulli  equation  along  a 

gne  de  courant  a  etc  presentee.  Dans  l’analyse  ma-  streamline  has  been  studied.  In  the  niatheiuatical 

theniatique  du  probleme,  l’oquation  de  Bernoulli  a  analysis  of  the  problem,  the  Bernoulli  equation  was 

ete  dei'ivee  a  partir  des  equations  fondamentales  du  derived  from  the  basic  equations  of  motion,  and 

mouvement  et  ensuite  siniplifiee  pour  les  quatre  cas  then  simplified  for  the  four  eases  under  consider* 

consideres  dans  cettc  etude.  Cependant  il  est  a  noter  ation  in  this  study.  However,  no  information 

qu’aueun  renseignement  concernant  I’energie  cin<5-  regarding  the  kinetic  energy  of  the  turbulent  flow 

tiijue  de  la  turbulence  de  menie  que  la  vraie  valeur  can  be  gained  by  this  use  of  the  Bernoulli  equation, 

<Ie  1’encrgie  perdue  localement  par  dissipation  vis-  nor  is  a  true  measure  obtained  of  the  energy  lost 

queusc  lie  peut  etre  obtenu  par  cette  equation  de  locally  through  viscous  dissipation.  The  results, 

Bernoulli.  Les  resultats,  montrant  les  valeurs  des  showing  the  values  of  the  pertinent  Bernoulli 

differents  tenues  de  l’equation  de  Bernoulli,  ont  6te  terms,  have  been  represented  in  graphical  form 

representes  sous  forme  de  graphique  pour  differen-  for  different  sections  along  a  streamline.  These 

tes  sections  le  long  d’une  ligne  de  courant.  II  ap-  results  lead  to  the  following  conclusions: 

parait  de  par  ces  graphiques  que  : 


La  plus  grande  partie  de  l’energie  perdue  a  lieu 
dans  la  region  du  iourbillon  ct  est  due  a  la  produc¬ 
tion  de  turbulence. 

»*  En  eonsiderant  uniquenient  l'ecou lenient  princi- 
*  pal,  il  apparait  quo  lc  terme  representant  le  trans- 
fert  d'energie  est  superieur  a  celui  representant  la 
f;  production  de  turbulence  niontrant  ainsi  que  la 
plus  grande  partie  de  l’dnergie  est  transferee  au 
iourbillon  et  faisant  ressortir  par  la  liienie  le  role 
joue  par  le  tourbillon  dans  les  phenonienes  de  sepa¬ 
ration. 

La  valeur  prise  par  la  sonime  des  terines  de  Ber¬ 
noulli  decroit  partout  a  l’exception  de  la  zone  du 
tourbillon  uu  voisinage  de  recouiemenl  principal  ou 
un  gain  local  d’energie  par  transfert  excede  la  va¬ 
lour  prise  pour  la  production  de  turbulence. 

L’upplicalion  du  tlieoreme  de  Bernoulli  a  un 
eeoulement  turbulent  est  certainement  une  verifica¬ 
tion  etilcace  d’une  analyse  globale;  elle  est  un  bon 
moyen  par  lequel  le  role  joue  par  le  tourbillon  peut 
ctre  eelairci  et  enfin,  elle  nous  fournit  sous  cer- 
taines  hypotheses  une  methode  de  determination  de 
ce  que  les  tensions  tangontielles  de  turbulence  de- 
vraient  etre. 


1.  Most  of  the  loss  of  energy  takes  place  in  the  re¬ 
gion  of  the  eddy  and  is  due  in  both  primary  and  sec¬ 
ondary  flow  to  the  production  of  turbulence. 

2.  By  considering  only  the  primary  flow,  it  is  seen 
that  the  transfer  term  is  higher  than  the  turbulence 
production  term,  thus  showing  that  most  of  the  energy 
lost  from  the  primary  flow  has  been  transferred  to  the 
eddy.  This  large  amount  of  energy  infused  into  the 
eddy  shows  clearly  its  role  in  the  separation  phe¬ 
nomenon. 

3.  The  usual  Bernoulli  sum,  in  general,  decreases 
everywhere  except  in  that  part  of  the  eddy  dose  to 
the  primary  flow,  where  the  local  gain  of  energy 
through  transfer  is  in  excess  of  the  local  turbulence 
production. 

4.  This  application  of  the  Bernoulli  theorem  for 
turbulent  flow  in  a  region  of  separation  provides,  un¬ 
der  certain  assumptions,  a  method  of  determining 
what  the  turbulence  shear  stresses  could  be;  it  is  cer¬ 
tainly  an  efficient  check  for  an  over-all  analysis,  and 
it  is  a  good  means  by  which  the  role  of  the  eddy  can 
be  elucidated. 


Repartition  <le-  Diets  llquides,  pour  une  parol  normnlc  /T/  Pattern  of  i stream  lines  for  normal  wall 
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Repartition  des  valcurs  mesurdes  /9 
de  l'elfort  de  cisalllement  turbulent, 
pour  un  dlsque  circuluire, 
dans  les  sections  correspomlant  U  1  <  l/Do  <  6 
Distribution  of  measured  turbulence  shear  stress 
for  circular  disli  in  sections  l  ^x/li •  <!) 


Oo 


10/  Repartition  des  valeurs  recherchdes 
de  l’effort  dc  cisaillement  turbulent, 
pour  un  disquc  circulalrc, 

dans  les  sections  correspondant  i  K  -r/Do  <  6 
Distribution  of  required  turbulence  shear  stress 
for  circular  disk  in  sections  1  ^  x/D,  ^  (I 
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